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Results are obtained for the decay rate of ortho and para states and for the hyperfine structure of 
the dimuonic system p,~ (dimuonium). We calculate next-to-leading order radiative corrections. 

The decay rate is strongly influenced by the electronic vacuum polarization in the far time-like 
asymptotic region and thus allows for a test of QED in a previously unexplored kinematic regime. 
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The production of exotic atoms and their properties have been discussed in a number of recent publications (see e.g. 
1^-^). Pionium, the bound system of a positive and a negative pion, has been observed in a beam target experiment 
with proton projectiles at the Serpukhov accelerator |^. Atomic bound states were also observed in decay modes of 
certain particles, e.g. positronium can be formed in the -s- positronium + 7 decay of the neutral pion [Q, and the 
TT p, atom can be produced in the reaction —> (tt p, atom) + :/ Q . 

The bound system consisting of two muons (dimuonium) can be produced in heavy nuclei inelastic scattering 
at high energies and in particle decays. The decay of the neutral meson into dimuonium has been investigated 
theoretically by L. Nemenov Q ( 77 ° —> dimuonium + 7 ). The formation of dimuonium in pion-proton collisions 
(7r“ + p —> dimuonium + n) and by photons on nuclei (j + Z dimuonium + Z) has been discussed by S. Bilenkii 
et al. 0 - For the direct production of dimuonium in muon-antimuon collisions, experimental difficulties associated 
with slow muon beams would have to be overcome. Another possible pathway for the production of the system, which 
we do not discuss in any further detail here, is the e“'"e“-annihilation (near or above the p threshold). Dimuonium, 
once produced, undergoes atomic decay (into energetically lower atomic states) and annihilation decay (into electrons 
and photons). Because the annihilation products are hard photons and relativistic electron-positron pairs, the decays 
could be investigated experimentally by established methods of particle physics. 

The lifetimes of low-lying S states in the dimuonic system lie in the 10“^^ s range. The decay products are hard 
photon pairs and relativistic electron-positron pairs, which can be easily detected. Because decay rate measurements 
can usually be performed with only a small number of events available, it is natural to begin the study of exotic 
systems with a detailed analysis of the decay channels and rates. 

In this Letter we obtain results for the decay rates of the n^Si ortho and the n^So para states of the dimuonic 
system for n = 1 and n = 2. We evaluate radiative corrections in next-to-leading order. We analyze the hyperfine 
structure of the system and obtain results for n = 1 and n = 2 states. 

The name “dimuonium” was proposed by J. Malenfant for the bound system of two muons in his pioneering 
investigations on the system ||^,^. In analogy to positronium we speak of ortho- and paradimuonium. We observe 
that a precise measurement of the orthodimuonium decay rate would allow for sensitive tests of QED. We emphasize 
that both the hyperfine structure and the decay rate of dimuonium are influenced by the electronic vacuum polarization 
in the far time-like asymptotic region (due to virtual photon annihilation processes). This regime does not yield any 
contribution in any other muonic or electronic bound system (see e.g. [ pd| , |12|). The electronic vacuum polarization 
contributes a relative correction of 1.6% to the decay rate of orthodimuonium, and thus would allow a test of QED 
in a previously unexplored kinematic regime. 
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The hadronic contribution to the vacuum polarization (for time-like below the pion threshold) can be found 
from a precise measurement of the orthodimuonium decay rate or the hyperfine structure. The correction enters at 
the level of one permille. 

We study here radiative corrections to the orthodimuonium and paradimuonium decay (Fig. ^jand respectively). 
We also investigate the hyperfine splitting (hfs) of the ground state (Fig. ^ and Fig. ||). Some of the contributions to 
the hyperfine splitting are closely related to corrections for the orthodimuonium decay, and certain corrections can be 
directly found in analogy with the positronium hyperfine splitting (Fig. |^), which has been considered e.g. in |p^ , p^ . 
The leading term for the hyperfine splitting is the Fermi energy Ep = (7/12) m, where m = rrie for positronium and 
m = rrifi for dimuonium (fi = c = 1). We list here the next-to-leading order contributions to hfs in the order in which 
they are depicted in Fig. |^. They are due to the vertex correction (fermion line), recoil corrections, vertex corrections 
in the annihilation part, vacuum polarization due to a loop of the constituent fermion (e“'"e“ for positronium, 
for dimuonium), and virtual two-photon-annihilation. These contributions, respectively, are given by (see IMD 
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For dimuonium, the small length scale of the system enhances the effects of the vacuum polarization corrections. 
The electronic vacuum polarization generates corrections to the wave function which enter at the order of a/n times 
the leading term as well as a large contribution due to an electron-positron vacuum polarization loop in the virtual 
annihilation diagram. The additional corrections specific to dimuonium, which are not present in positronium, are 
depicted in Fig. ^ The contributions are given by the (electronic) vacuum polarization corrections to the wave 
function in the transverse photon exchange diagram, by a vacuum polarization insertion in the transverse photon line, 
by vacuum polarization corrections to the wave function in the annihilation diagram, and by electronic and hadronic 
vacuum polarization insertions in the virtual photon line (annihilation diagram). Respectively, we obtain the following 
results for the additional contributions specific to dimuonium, 
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It should be noted that the hrst three contributions in Eq. (||) are state dependent, i.e. they do not exhibit the 
familiar 1/n^ scaling. Because atomic momenta are of the order it is sufficient to carry out the calculations 

with non-relativistic wave functions (cf. |l^ ). The numerical methods employed in [p^ for the treatment of the non- 
relativistic Green’s function are also used for the calculations presented in this Letter. The results have been checked 
by independent evaluations using the computer algebra system Mathematica ||1^, and will be described further in 


In the discussion of the hadronic vacuum polarization, we follow an approach outlined in fl^ ] . We thereby estimate 
the uncertainty due to our model as 11 % of the total hadronic vacuum polarization contribution. This is sufhcient 
because the hadronic contribution is small compared to other corrections, and because higher order QED radiative 
corrections are expected to contribute in the same order as the uncertainty of our model of hadronic vacuum polariza¬ 
tion. The main contribution to the hadronic vacuum polarization originates from the p meson pole in the pion form 
factor. 

The real part of the sum of all the next-to-leading order corrections for dimuonium is 
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For the 2S state, the contributions sum up to £’^fg(2S') = (o/tt) 0.556 • Epji. We obtain the following theoretical 
values for dimuonium hyperfine splitting, Ehfs(lS') = 4.23283(35) • 10^MHz and Ehfs(2S') = 5.28941(34) • 10® MHz, 
where we estimate higher order quantum electrodynamic corrections to enter at the level of 5 % of the result obtained 
for next-to-leading order corrections. The hyperfine structure of dimuonium could be measured in an infrared laser 
field which mixes the two states and thus modihes the decay rate of the statistical sample. 

The main decay channel for orthodimuonium is given by one-photon annihilation into an electron-positron pair. 
The first-order result can be extracted from the imaginary part of the electronic vacuum polarization contribution 
(annihilation diagram) to the hyperhne splitting (Eq. (^)). The result is F^®)(n®S'i) = a®m^/(6n®). Corrections to 
this result (see Fig. ||) originate from similar diagrams as those for the hyperfine splitting. Some of the calculations 
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are therefore closely related to those for the hyperfine splitting. The next-to-leading order corrections are due to 
radiative corrections to /r lines, due to the muonic, electronic and hadronic vacuum polarization, due to the wave 
function correction caused by electronic vacuum polarization, due to radiative corrections in the electron-positron- 
line and bremsstrahlung corrections, and three-photon decay. We obtain for the respective contributions, 
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For the 2S state, we obtain the result Ar(2^S'i) = {a/ir) 3.60(4) ■ r(°^(2^S'i). 

Corrections to the orthodimuonium decay rate have been investigated previously In |^, only the decay of 

orthodimuonium into a single electron-positron pair and into an electron-positron pair and a photon is taken into 
consideration. Three-photon decay, however, constitutes an important correction to the overall decay rate (last term 
in Eq. (|)). 

There is also a numerical discrepancy between the result reported in and our result for the contribution of 
the electronic vacuum polarization correction to the wave function (fifth diagram in Fig. fifth term in Eq. (^). 
We presume this disagreement can be traced to the fact that J. Malenfant has calculated this correction with free 
Green functions, whereas the evaluation in this work is done using bound Green functions. We hold the view that it is 
indispensable to use bound Green functions for the respective correction, because the atomic momentum in dimuonium 
is of the order of a to^/ 2 « 0.75 TOe. This momentum is close to the mass of the loop particles (electrons and positrons) 
of electronic vacuum polarization. These particles determine the radius of the Uehling potential. Therefore, the decisive 
region of the momentum integration for this correction is about TOe, where the bound Goulomb Green function cannot 
be approximated by the free Green function (because the effect of the binding Coulomb potential, in momentum 
space, is inversely proportional to the square of the momentum transfer). 

Our final theoretical results for the decay rates of the low lying ortho states are r(1^5'i) = 1.79073(77) • 10“^^ s and 
t(2^S'i) = 1.43305(59) • 10“^^ s, where we have assumed that higher order corrections contribute at the level of 5% 
compared to next-to-leading order corrections. 

In this Letter, we also investigate the decay modes of the para system. The main decay channel of paradimuonium 
is given by the annihilation into two (real) photons. The first-order result for the decay rate can be extracted from the 
imaginary part of the next-to-leading order correction to positronium hyperfine splitting (Eq. (Q)). The imaginary 
part in Eq. (|lj) is entirely due to the virtual two-photon annihilation of the para state. The decay rate in leadiM 
order can be inferred as T^^\n^So) = a® to^/( 2 n^). The next-to-leading order corrections are depicted in Fig. |j. 
They are given by radiative corrections in the muon lines, by the wave function correction due to electronic vacuum 
polarization and by two-photon annihilation processes with production of an additional electron-positron pair. We 
obtain for the respective contributions, 
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For the next-to-leading order corrections to the decay of the 2^5'o state, we obtain Ar(2^5'o) = {a/n) 4.65 -F^^^ (2^5'o). 
The theoretical values for the decay rates of para states are r(l^S'o) = 5.9547(33) • 10“^^ s and t(2^5'o) = 4.7653(25) ■ 
10“^^ s, where we have estimated the magnitude of higher order QED corrections as 5 % of the corrections considered 
above. 

The last two diagrams in Fig. which correspond to the decay paradimuonium ee + j, contribute at the 1 % level 
to the decay rate of the para system. The respective correction contains the large logarithm of the muon to electron 
mass ratio (the 4/3 ln(2TO^/TOe) — 16/9 term in Eq. (^). This contribution has been evaluated in two independent 
ways. One of the methods employed is standard S-matrix formalim, using the computer algebra system Mathematica 
0 for the analytical calculations. For a simplified treatment, we observe that the decay rate of paradimuonium into 
two photons can be written as the imaginary part of the virtual two-photon annihilation (last two diagrams in Fig. ^). 
In order to obtain a result for the paradimuonium —> ee -I- 7 decay, we consider vacuum polarization insertions in the 
photon lines. The result can then be obtained by considering correction to the photon propagator due to electronic 
vacuum polarization, and by analyzing the decay of paradimuonium into two photons, one with vanishing mass, the 
other having a mass that corresponds to the sum of the four-momenta of the emerging electron-positron pair. More 
details of this calculation will be presented elsewhere ]T^ . 

In conclusion we obtain the value of Ehfs(lS') = 4.23283(35) • 10^ MHz for the hyperfine splitting of the l^ states 
in dimuonium. This frequency is lying in the infrared region and can be produced with e.g. diode lasers. For the 
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lifetime of the ortho state we obtain r(l^S'i) = 1.79073(77) ■ 10“^^s. The para state has a lifetime of 
t(I^S'o) = 5.9547(33) • 10“^^ s. Some of the corrections to the decay rate are state dependent and break the 1/n^ 
scaling of the corrections found in the positronium system. Because they can be expressed as a correction to the 
wave function at the origin, caused by the Uehling potential (e.g. the fifth diagram in Fig. these corrections are 
more sensitive to details of the vacuum polarization potential around the origin than integral values like the energy 
shift. Due to the small length scale of the muonic system, hadronic vacuum polarization corrections also contribute 
at the level of one permille to the decay rate of the ortho state. It should be noted that the enhancement of the 
orthodimuonium decay rate by 1.6 % due to the electronic vacuum polarization in the far time-like asymptotic region 
allows for a sensitive test of QED (the contribution is given by the third term in Eq. (^). This kinematic region 
cannot be explored in other atomic systems. 
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FIG. 1. Known corrections to the hyperfine splitting of S states of a dileptonic system. Springy lines denote transverse 
photons. 
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FIG. 2. Additional corrections to the hyperfine splitting specific to dimuonium. Bold fermionic lines denote muons, thin 
lines denote electrons. Dashed lines denote Goulomb photons. 


5 









-L.O. 







+ symmetrization 


FIG. 3. Feynman diagrams for next-to-leading order corrections to the decay rate of orthodimuonium. The subtraction of 
lower order (L.O.) terms is necessary for some contributions in order to prevent double counting. 
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FIG. 4. Feynman diagrams for next-to-leading order corrections to the decay rate of paradimuonium. 
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